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Abstract 

In this paper we present an investigation of numerical Monte Carlo sim- 
ulations of the diffusive shock acceleration in the test particle limit. Very 
high gamma flow astrophysical plasmas, have been used, from j up ~ 50 up 
to "f U p ~ 1000 , which could be relevant to the suggested models of AGNs 
Jets and their Central Engines as well as the ultra-relativistic shock particle 
acceleration in Gamma Ray Burst (GRB) fireballs. Particularly the energy 
gain per shock crossing and the time constant for the above high relativis- 
tic processes is numerically calculated. We explicitly find a considerable 7 2 
energy boosting in the first shock cycle, and in all subsequent shock cycles 
the particle energy is multiplied by a large factor. Also a noted accelera- 
tion speed-up for the same acceleration process has been observed. Both 
of those results are connected with theoretical suggestions that a significant 
enhancement of the acceleration process is possible due to relativistic effects 
not present at lower plasma flow speeds. The acceleration speed-up and the 
mean energy-gain per shock cycle found could efficiently justify the theory 
of the origin of Ultra High Energy Cosmic Rays (UHECR) from the sites of 
GRBs. 

1 Introduction 

Ultra-relativistic flow velocities are found in many astrophysical objects such as in 
Gamma Ray Bursts (GRB) as well as in AGN jets and their central engines where 
accretion is observed. The current understanding of multiwavelenght GRB after- 
glow observations, indicates that they are produced at the relativistic shocks, formed 
shortly after a violent explosion takes place where a blast wave impacts with the 
surrounding interstellar medium. Consequently shock waves form (Meszaros & Rees, 
1993), and eventually relativistic diffusive shock acceleration takes place. Theoretical 
works of Vietri (1995,1997 and 1998) and Waxman (1995) suggested that diffusive 
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acceleration at non-relativistic shocks is not fast enough to produce the highest en- 
ergy particles observed, but they considered that a potential relativistic supesonic 
environment could succeed to do so. One of the first computational works (Quenby 
& Lieu, 1989)) concerning relativistic (7 ~3) diffusive acceleration with application 
to AGN relativistic flows reported a considerable enhancement in the acceleration 
rate of a factor~13. In this present work we present a numerical investigation for 
diffusive acceleration in ultra-relativistic shocks with application to GRBs theoretical 
models. 

The standard analytical theory of diffusive shock acceleration operating in one 
dimensional astrophysical non-relativistic flows, where u±, u 2 are the non-relativistic 
upstream and downstream flow speeds and the particle's velocity is almost equal to 
c, yields a differential spectrum dn/dp oc p~ a , where a = r + 2/r + 1 and r = U\ju 2 
and the analytically derived acceleration time constant is given by, 



Ui — U 2 Ui u 2 

where Ai, A2 are the particle's mean free paths, upstream and downstream respec- 
tively. It is known that the fractional energy increase for a single shock crossing, for 
relativistic flows, is ~ (1 — V 2 / c 2 )~ a5 ) where V is the relativistic difference {ui — u 2 ) 
in flow velocities across the shock. Considering again relativistic flows and diffusion 
theory calculations for the cycle-time constant, it is found that the latter does not 
change from the non-relativistic case, and so we could expect a speed-up of acceler- 
ation, relative to the non-relativistic estimate. Because relativistic flows create large 
anisotropics in the plasma flow a detailed computational approach is necessary in this 
high 7 regime. Several studies have been made by many authors, but still, there is 
a controversy about the rate of particles energy gain and the time constant of the 
acceleration. These controversies need a solution as they are crucial for models of 
AGN jet cosmic ray acceleration and UHECR production and acceleration in GRBs 
(eg Vietri). In this present work we present a Monte Carlo simulation, investigating 
the diffusive acceleration in ultra-relativistic parallel shocks for values up to 7 ~ 990, 
with application to GRBs relativistic flows. We find a considerable acceleration rate 
enhancement and our results show explicitly that the energy boosting, scales as 7 2 in 
the first shock passaging cycle. After that, the particle energy is multiplied by a large 
number, gaining large amounts of energy which could still account for the theory of 
UHECR. These results support theoretical predictions (eg. Vietri 1995) concerning 
UHECR origin and acceleration from GRBs' highly relativistic plasma flows. 
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2 Numerical method 



In this Monte Carlo code we consider both isotropic large angle and pitch angle 
diffusion which is calculated in the respective plasma rest frames, where particles are 
allowed to scatter in the respective fluid frames towards the shock. We consider only 
parallel shock (0 = 0, where 9 is the angle between the shock normal and the magnetic 
field), either because that is the field configuration, or turbulence removes "reflection" 
at the interface. A guiding centre approximation is used where the particle trajectory 
is followed in one-dimensional space along the x axis. A relativistic transformation is 
performed to the local plasma frames each time the particle scatters across the shock 
following it according to particle jump conditions and it is made leave the system 
from the moment that it 'escapes' far downstream at the spatial boundary or if it 
reaches a well defined maximum energy. The particles (~10 7 ) of weight equal to one, 
are injected far upstream at a constant energy of high gamma, which supposes that a 
pre-acceleration of the particles has already taken place. They left to move towards 
the shock where along the way they collide with the pressumed scattering centers 
and consequently as they keep scattering between the upstream and downstream 
regions they gain each time an amount of energy. The compresion ratio is allowed 
to have the value of 4, for immediate comparison with the non-relativistic values. 
A main characteristic of the code is that a particle splitting technique is used, in 
order for the statistics to be as efficient as possible. The main notion of particle 
splitting is the fact that one can consider that each particle represents instead a 
large number of particles. When a number of particles escapes through the defined 
spatial or momentum boundaries during the acceleration process, we replace these lost 
particles with new ones -in order to have again the same initial number of particles- 
but in such a way that as their weight is decreased, the number of particles remains 
almost constant throughout the simulation. The probability that a particle will move 
a distance z along the field lines at pitch angle 9 before it scatters is given by the 
following expression, 



For the case of pitch angle diffusion scattering, the new pitch angle 9' is calculated 
by the trigonometric formula, 



Prob(z) ~ exp{— z/\\cos9\) 



(2) 



cos9' = cos9Vl — sin 2 59 + sin59Vl 



— cos 2 9coscj> 



(3) 
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Figure 1: Left-The ratio of the computational to the non-relativistic analytical time con- 
stants for large angle scattering versus the Lorentz ^ up Qow.Right-The equivalent plot for 
pitch angle diffusion. As it is seen, for both plots there is a maximum substantial 'speed 
up' of a factor ~5 and ~200, respectively. 

where <f>e(0,2n) is the azimuth angle with respect to the original momentum direction. 
We note here that for the case of highly relativistic flows, the definition of particles 
pitch angle diffusion is scattering within an angle ~ I/71, where 7 X is the upstream 
gamma measured in the shock frame (Gallant & Achterberg, 1999). This important 
condition is explicitly included in our Monte Carlo simulations during the relativistic 
particle shock acceleration process, having as a main aim to investigate and compare 
our results with controversial theoretical predictions which depend whether either 
large angle or pitch angle diffusion operates in highly relativistic flows. 

3 Results 

The aim of this work is to investigate the energy boosting of particles accelerated in 
higly relativistic shock waves in comparison to the shock cycle and to study the cycle 
time constant for the particle acceleration, by simulating highly relativistic shock 
formations which potentially occur in GRBs. The results shown in the figures, have 
been recorded in the downstream side, in the shock frame and concern large angle 
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Figure 2: Left: For large angle scattering we calculate the mean energy of the particle versus the no. 
of shock crossings, employing an upstream flow gamma factor 7=50. We may observe the efficiency 
of the 7 2 boosting in the first cycle and the large energy multiplication in all subsequent crossings. 
Right: For pitch angle diffusion the particle energy gain appears (normally) to be less high. 



large angle scattering 



pitch angle diffusion 




no. of crossings 



4 6 
no. of crossings 



Figure 3: Left: In this plot we use large angle diffusion and the upstream flow has the value of 7=200. 
It is seen again the efficiency of the mean energy gain in every shock crossing Right: For pitch angle 
diffusion, the mean energy gain per cycle for 7=200 versus the shock crossing, gives a considerable 
energy gain efficiency mostly in the first shock cycle crossing (7 2 ). 



5 



large angle scattering 



jamma = 990 



2 3 4 

no. of crossings 



Figure 4: In this plot we use 7=990 for large angle diffusion. One may observe that the 
energy boosting efficiency is explicitly present in the first shock cycle. 
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Figure 5: For pitch angle diffusion, the mean energy gain per cycle for 7=990, to be 
compared with the square of the upstream flow gamma factor 7 2 for the first shock cycle. 
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isotropic scattering and pitch angle diffusion for a various number of upstream Lorentz 
7 values. Figure 1, shows for both large angle scattering and pitch angle diffusion, the 
ratio of the computational time constant to the non-relativistic analytical acceleration 
time constant, as a function of the upstream Lorentz factor 7 np plasma flows. We note 
here that in the case of the pitch angle diffusion, we allow for the mean free path of 
particle to scatter ~ ir/2 for the calculation of the analytical time-expression. We can 
observe that substantial 'speed-up' of a factor ~ 5 for large angle scattering is given 
and a considerable 'speed-up' of a factor ~ 200 is noted for pitch angle diffusion. The 
figures 2-5 show the mean energy gain of the accelerated particles versus the shock 
crossing where a substantial 7 2 boosting is found, explicitly for the first cycle. After 
that, and for the following shock crossings, the energy of the particle is multiplied 
by a large factor. Our first indication ( / ~f 2 ) is in good aggreement with Gallant & 
Achterberg (1999) predictions who found 'that the energy of the particle for the first 
shock cycle is increased by a factor of 7 2 p '. On the other hand though, they found 
that ...'in all subsequent shock crossing cycles the particle energy doubles' and that 
due to the fact that the particle does not have time to re-isotropise upstream before 
being overtaken by the shock. Our results show that the particle energy is multiplied 
by a large factor. This could be justified from the moment we could assume that 
the acceleration speed-up found, gives enough time for the particle to re-isotropise 
upstream before being overtaken by the shock. 

4 Conclusions 

A Monte Carlo numerical investigation has been reported in the test particle limit 
of parallel diffusive shock acceleration. Very high gamma flow astrophysical plasmas 
have been used, from ^ uv ~ 50 up to 7„ p ~ 1000 which are relevant to highly relativis- 
tic astrophysical supersonic plasma flows in GRBs. A dramatic 7 2 energy boost in the 
first shock crossing cycle has been observed and a considerable acceleration speed- 
up has been found for both isotropic large angle and pitch angle diffusion. Both 
these results support Vietri and Waxman GRB models and theoretical predictions 
for the model of Ultra High Energy Cosmic Ray origin from GRBs. The surprising 
implication of those results is that a modeler could consider that, indeed there is a 
distinctive case showing that a 7 2 energy gain versus first shock crossing and a large 
energy multiplication in all other subsequent cycles is indeed happening in highly 
relativistic shock waves, 'making space' to the particles in order to gain rather very 
large amounts of energy in few shock cycles. On the other hand the above impli- 
cations, along with our numerical calculations, finding a considerable acceleration 
speed-up could sufficiently account for the theoretical predictions implying UHECR 
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origin from GRB relativistic shock configurations. However these crucial issues, need 
further detailed investigation as many other parameters need to be included within 
the simulation codes such as shock obliquity, non-linear processes and energy losses in 
order for potential models to be as realistic as possible with the actual relativistic flow 
of GRBs. There is already under way a similar work concerning the above parameters. 
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